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ABSTRACT 


This  document  outlines  a  guidance  equation  mechanization 
for  use  with  ballistic  missiles.  This  technique  does  not 
require  pre-launch  targeting  computations  and  is  therefore 
especially  suitable  for  mobile  weapons  systems. 
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SECTION  I 


INTRODUCTION 


The  development  of  a  guidance  technique  suitable  for  use  with  any  ballistic 
missile  is  presented  in  this  document.  This  technique  is  particularly  suita¬ 
ble  for  weapon  systems  requiring  maximum  flexibility  (mobile  systems  re¬ 
quiring  a  minimum  operational  targeting  time).  While  an  example  of  the 
basic  equations  required  for  mechanization  is  given^  relatively  routine  as¬ 
pects  of  the  mechanization  (such  as  data  reduction  and  first  stage  and  vernier 
guidance)  are  not  discussed.  This  mechanization  is  obviously  more  compli¬ 
cated  than  that  which  would  be  required  if  the  equations  were  not  self  con¬ 
tained  (i.e. ,  if  targeting  operations  are  carried  out  at  the  launch  site). 

The  equations  of  motion  which  have  been  used  for  the  explicit  solution  on  the 
free  flight  trajectory  are  derived  from  the  assumptions  that  the  earth's 
gravitational  field  varies  with  the  inverse  square  law  and  that  there  is  no 
atmosphere  (see  Section  3),  The  inadequacy  of  this  assumption  can  be  com¬ 
pensated  for  by  empirically  determined  target  offsets,  but  a  considerably 
superior  method  which  does  not  require  any  targeting  operations  at  the  op¬ 
erational  site  is  developed  in  Section  4. 

The  only  inputs  required  by  the  guidance  equations  are  latitude,  longitude, 
and  altitude  (above  sea  level)  of  the  launcher  and  target.  The  guidance  equa¬ 
tions  then  compute  all  other  quantities  required,  including  launch  azimuth. 
The  guidance  inputs  need  to  be  specified  sufficiently  prior  to  launch  to  permit 
the  launch  azimuth  to  be  implemented  (computation  time  is  negligible).  With 
the  possible  exception  of  hardware  requirements  there  is  no  other  reason  for 
a  delay  between  the  receipt  of  information  regarding  target  and  launcher  lo¬ 
cations  and  launch  of  the  vehicle. 


Simulation  has  shown  that  tha  on#  r  miss  distaacs  da#  to  ths  inaccuracias  in 
this  technique  can  easily  be  reduced  to  200  feat.  Considerably  batter  accu* 
racy  may  be  obtained  if  empirical  coefficients  are  determined  carefully. 
The  above  error  does  not  include  the  effects  of  gravitational  anomalies » 
hardware  inaccuracies  or  measurement  errors. 


} 
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SECTION  2 


DEFINITIONS 


The  following  definitions  are  given  to  aid  in  understanding  the  discussion 
which  follows: 

Explicit  Guidance  is  a  generic  term  for  the  class  of  guidance 
equations  which  are  a  direct  soluticm  of  equations  of  motion. 

The  equations  of  motion  which  are  used  may*  or  may  not»  rep¬ 
resent  reality  with  sufficient  accuracy  to  satisfy  mission  objec¬ 
tive;  there  is,  however,  a  restriction  in  the  latter  case  that  the 
inaccuracy  (which  is  compensated  for  by  empirical  or  semi- 
empirical  methods)  must  be  extremely  small. 

Simulation  is  a  computer  (usually  high-speed  digital)  program 
whlcti^s  M  accurate  representation  of  the  earth,  atmosphere, 
missile  engine  and  autopilot,  etc.,  as  well  as  the  guidance 
equations.  This  simulation  then  "flies"  the  missile  and  deter¬ 
mines  performance  capability,  workability  of  the  guidance 
equations,  etc. 

Targeting  is  the  utilization  of  a  simulation  to  define  and  verify 
any  empir i cal  constants  that  may  be  required  by  the  guidance 
equations.  For  many  problems  efficient  use  of  the  simulation 
to  obtain  the  empirical  constants  requires  the  use  auxiliary 
computer  programs;  this  kind  of  effort  is,  of  course,  part  of 
targeting.  The  word  targeting  is  sometimes  applied  to  opera¬ 
tions  (carried  out  at  the  operational  site)  which  utilize  the  empir¬ 
ical  constants  obtadned  by  the  process  defined  above  as  targeting. 

Basic  Inertial  Coordinates  ^  are  earth  centered  with  4 

passing  through  ^e  nor&  pole.  These  coordinates  and  their 
derivatives  are  scaled  so  that  the  gravitational  constant  niulti- 
plied  by  the  mass  of  the  earth  (GM)  is  unity.  At  some  time  dur¬ 
ing  the  pre -launch  operation,  ^e  gyros  are  uncaged.  The  inertial 
coordinate  system  is  defined  at  this  instant  and  the  quantity  T« 
represents  the  time  interval  which  has  elapsed  since  this  instsnt. 
These  coordinates  would  not  necessarily  be  used  in  an  actual 
mechanization,  particularly  for  an  inertial  system. 
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SECTION  3 

A  KEPLERIAN  SOLUTION 


Coordinate  Systems 


It  is  necessary  to  relate  the  position  and  velocity  of  the  missile  and  the 
target  by  means  of  an  inertial  coordinate  frame.  Once  this  has  been 
done,  any  appropriate  criteria  can  be  used  to  develop  expressions  for 
pitch  and  yaw  velocities.  Criteria  which  have  been  used  in  the  past  and 
have  proven  to  be  very  satisfactory  are  outlined  below. 


The  missile,  target  and  geocentroid  define  a  plane.  The  yaw  component 
of  velocity  is  the  component  of  missile  velocity  normal  to  this  plane 

and  the  pitch  component  of  velocity  is  the  component  of  missile  velocity 

• 

in  this  plane  and  parallel  to  the  radius  from  the  geocentroid  to  the 
missile.  Pitch  angles  are  measured  positive  up  from  the  normal  to 
this  radius.  (See  Figure  1  and  Appendix  A. ) 


The  subscript  D  is  used  to  denote  burnout  conditions.  Adequate  estimation 
of  vjijj  is  straightforward  (see  Appendix  A). 

Orbital  Parameters 

The  value  of  is  defined  as  a  compromise  between  minimum  re-entry 
dispersions,  re-entry  heating  problems,  and  performance.  Actually  the 
numerical  coefficients  in  the  definition  of  F^  (Appendix  A)  can  be  chosen 
to  provide  a  good  approximation  to  any  reasonable  functional  relationship 
for  F^.  The  desirability  of  steering  to  a  given  F^  is  related  to  re-entry 
effects  and  is  discussed  in  Section  4. 


I 
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The  expression  for  the  velocity  necessary  to  travel  from  any  point  in 
space  to  any  other  point  (assuming  no  atmosphere  and  an  inverse  square 
gravitational  field)  can  be  derived  from  Figure  Z  and  the  well  known 
orbit  equation 


Ar  a/  r  a 

where  a  is  the  semi-major  axis  of  the  orbit.  From  Figure  2 

(2a  -  =  [r^  sin  4^  -  (2a  -  r^)  sin 

+  [rj  -  (2a  -  rj)  cos  2rj  -  r^  cos  4^]^ 


or 


1  + 


r2(l  -  cos  4^) 

r  j(l  +  cos  2rjT  •  r^Tl  +  cos  (4^  +  ZEj  )) 


(2) 


Combining  Equations  1  and  Z 


‘  '  'iV  2“/  ■'l 


2(1  -  cos  4^) 


—  (1  +  cos  2r. )  -  cos  4^  -  cos  (4^  +  zr. ) 


(3) 


Defining 


6  H  rV^ 
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and  choosing  as  the  target  and  as  burnout 


6^  =  — 
T  Tr 


2(1  -  cos 


- —  (1  +  cos  2r^)  -  cos  «  cos  (44j^  + 


T  Vr 


from  Equation  1 


From  orbital  kinematics  (Reference  2) 


€  sin  E  =  >/6(2  -  6)  sin  r 


€  cos  E  =  6  -  1 


Tpp  =  a'^  ^  [E^  -  Ej  -  <  sin  ^2  ^ 


c  =  -  1)  +  6^(2  -  6^)  sm 


1  -  6^(2  -  6^)  cos 
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and 


6 


D 


=  2 


r 


D 


a 


V 


D 


An  expression  for  may  be  found  from  Equation  3  as  follows: 


2  *  D 

•g—  (1  -  cos  =  —(1  +  cos  ZFjj)  -  cos  -  cos  +  ZFjj) 


since 


cos  (4<jj  +  =  cos  »|<jj  cos  ZFjj  -  sin  sin  ZF^ 


Then 


sin  2Fp  =  -  Qj(l  +  cos  2Fjj) 


(11) 


where 


71  -  ‘•'D 

=  sin  .|<jj 

2(1  -  cos  41J5) 

02  = 


(U) 
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Squaring  Equation  1 1 


sin^  ZFjj  =  (Q^  -  Qj)'^  -  Z(Q^  -  Qj  )Qj  cos  cos^ 


or 


(l  +  Qj)cos^  ZFjj-  [2(Q2  -  Qi)Qj]  cos  2r  +  ((Q2  -  Qj  -  l)  =  0 


Then 


cos  zr 


D  ^ 


2(0^  -  Qj)Qj  ±  /4qJ(Q2  - 

q/  -M 

('  *  °‘\) 

(Q,  -  q/  -  1 

2| 

[l  .  Qf) 

0|(0^-0l)t'/|  -Q^t3Q,Q^ 

1  +  Q? 


Obviously  the  sign  of  the  radical  depends  on  whether  a  lofted  or  shallow 
trajectory  is  used  (the  radical  vanishes  at  ininimum  energy).  Since  a 
slightly  lofted  trajectory  will  be  used  (see  Section  4) 


cos  = 


Ql(Q,  - 


>/l 


1  -  q;  + 


1  +  Q. 


(13) 


A  convenient  specialized  solution. of  Kepler’s  equation  (Eq.  9)  for  ballistic 
missile  trajectories  can  be  developed  as  follows: 


€  cos  E  =  1  -  — 
a 
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Obviously  £  s  0  at  perigee,  and  since  a  is  less  than  the  radius  of  the  earth 
tr/2  <  E  <  3ir/2  throughout  the  trajectory.  However  if 


are  evaluated  in  the  first  quadrant  then  Equation  9  becomes 
Tpp  =  +  Ejj  +  €(iin  Ej^  +  sin  E^)J 


4 


(14) 


(15) 
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SECTION  4 


NON-KEPLERIAN  EFFECTS 


Introduction 

The  two-body  point  mass  central  force  field  problem  can  always  be  reduced 
to  an  equivalent  one-body  problem  (Reference  2).  For  an  inverse  square 
force  field  this  becomes  the  *'Keplerian  problem,  **  and  more  or  less  com¬ 
plete  analytic  treatments  have  been  known  for  centuries.  For  practical 
purposes  on  a  ballistic  missile  trajectory  the  solution  to  the  Keplerian 
problem  is  in  error  due  only  to: 

a.  Atmospheric  f  rces  at  re-entry 

b.  The  fact  that  the  earth's  gravitational  field  is  not  exactly 
inverse  square,  but  has  other  latitude  dependent  terms 

c.  Local  gravitational  anomalies,  and 

d.  Solar  and  lunar  gravity. 

Fortunately  these  effects  are  quite  small. 

Atmospheric  Effects  at  Re-entry 

For  practical  purposes  the  only  non-trivial  complications  arising  from 
inclusion  of  a  target  altitude  are  due  to  atmospheric  forces  at  re-entry. 

The  principal  result  of  these  atmospheric  forces  is  to: 

a.  Increase  the  time  required  to  reach  the  target  altitude 

b.  Cause  the  hit  point  to  be  less  far  "down  range"  in  the 
orbit  plane,  and 

c.  Impact  displacement  due  to  atmospheric  (earth) 
rotation. 

It  is  convenient  to  combine  (a)  and  (c)  to  obtain  an  "effective"  time  differ¬ 
ential  (which  is  quite  small).  Since  has  been  defined  as  a  function  of 
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(sec  Fig.  A-3)  the  two  important  re-entry  variables  (V  and  F)  have 
been  reduced  to  a  one  parameter  family.  With  this  simplification  it  is 
possible  to  compensate  for  re-entry  effects  (for  altitudes  of  up  to  one 
mile)  with  an  rms  error  of  about  100  ft  with  the  following  empirical 
equations: 

‘•'re  ""  ^  ^  ^  ^  S 

»  (kJ'*  ♦  >4^  *  k}Jvt)  (16» 

Equations  16  and  17  are  terms  in  ^nd  Tp  (see  Appendix  A). 

It  is  of  course  desirable  to  have  ^ud  as  small  as  possible  which 
would  require  a  very  lofted  trajectory.  However  when  this  effect  is 
balanced  against  the  requirements  of  heating  and  trajectory  sensitivity 
(miss  coefficients)  the  resulting  compromise  is  a  family  of  trajectories 
slightly  more  lofted  than  the  "minimum  energy"  family. 

Oblate  Potential 

The  earth's  potential  field  is  usually  represented  as 

GM„  it  ,  \  , 

V  =  -  1'^  +  J “  *i*>  6|  +  D  £(ain  9)  +  •  •  • 

o  r  '  r 

where  J  and  D  are  constants.  The  Keplerian  analysis  assumes  that  J,  D, 
and  all  higher  order  terms  vanish,  an  assumption  that  would  give  intoler¬ 
able  error  if  not  compensated  for.  The  basis  for  this  oblate  potential 
compensation  is  the  perturbation  development  by  Wheelon  in  Ref.  3. 
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The  D  term  is  about  one  tenth  percent  of  the  size  of  J  term  (not  one  percent 
as  stated  in  Ref.  3)  and  the  higher  order  terms  are  even  smaller,  and  since 
the  error  in  the  J  perturbation  is  larger  than  this  (due  to  neglect  of  the 
second  order  terms  in  the  perturbation  development)  primary  consideration 
will  be  given  to  the  first  order  perturbation  of  the  J  term.  All  residual 
errors  will  then  be  lumped  together  and  can  be  treated  empirically. 

Wheelon  develops  (Ref.  3Eqs.  65  -  69)  the  following  expression  for  cross 
range  miss  (in  present  notation); 


Wheelon  now  sets  r  constant  and  integrates.  A  more  general  expression 
is 


2Jr. 


6L  = 


I'd 


2 . .  sin  A^  cos  0^  [sin  0j^(l  -  cos 

+  C08  cos  1 


[Wheelon  has  ^2  vary 

slightly  depending  on  the  trajectory  used)  can  now  be  developed  empirically. 
For  ICBM  trajectories  of  tactical  interest  the  simple  expression 

f,  =  .  K,«l)  +  (Kj  +  K3.l'+*  '“O' 


is  very  satisfactory. 
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Since  6L  i«  effectively  a  change  in  the  deeired  plane 


\d  Tj  sin  (|»jj 


An  expression  similar  to  £q.  19  can  be  developed  from  Wheelon^s 
analysis: 

« =  ^--‘7'--°-^''  1‘3‘d  V  *  •  °12  ♦  V  j|  <«' 


6q  sin  cos 


where 


^  ^ 
+  cos^A^  cos^ 


The  functional  form  of  f^  is  chosen  to  be 


The  function  f^  arises  from  a  ^  term  and  thus  no  neat  semi^empirical 
relationship  can  be  found  by  a  constant  r  integration.  A  satisfactory 
representation  for  ICBM  trajectories  of  tactical  interest  has  been  foui^  to  be 


,,  .  (...^  .  c..^  co.^  A.)  ^0  *  kJ«/) 

2 

-  2  sin  coe  cos  +  kJ*  i|**  +  K*'  ♦*  J 


Equations  22,  23  and  24  can  be  combined  with  the  relationship 


J  = 


(25) 


To  obtain  a  correction  to  (see  Fig.  A-2). 

The  oblate  potential  also  changes  the  time  of  free  flight  from  that  given  by 
the  Keplerian  solution.  This  variation  can  be  approximated  empirically  by 
the  expression 


Tj  -  (K^  .  K,  e  J  t  (K^  t  Kj  »l)+*  .  k, 

*2 

-  il'  (K^  + 

Local  Gravitational  Anomalies 

Local  gravitational  anomalies  perturb  the  trajectory  not  only  in  free  flight 
but  in  powered  flight  as  well.  Unfortunately  these  anomalies  have  not  been 
mapped  completely  enough  to  allow  a  definitive  analysis  of  their  effect  on 
trajectories  of  general  interest.  Once  this  mapping  becomes  available  there  , 
is  no  conceptual  reason  why  these  effects  could  not  be  included  in  an  empiric 
cal  way,  although  the  amount  of  information  that  would  have  to  be  stored 
might  be  prohibitive . 


» 


Solmr  and  JLm»r  Gravity 

Solar  and  lunar  gravitational  fields  perturb  the  trajectory  whenever  the 
vehicle  and  the  geocentroid  are  not  equidistant  from  these  bodies.  The 
perturbing  acceleration  is 


where 


GM 

pert--^ 


GM  GM  /  ^  GM  2  6R 

If 


6R  =  R-  -  R,  <  radius  of  earth 
Z  1 


For  the  sun 

4 


^pert 


sun 


2«^  6R 
e 


where 


*  l°/day  «  2(10"^)  -i— 
6R«  2(10^)  ft 


■o 


pert 


»  1.6(10"Sft/«ec^ 


sun 


( 
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Since 


S  <ia  xi  =  0.8(10"S(2000)^*  3  feet 

pert^un"^  pert  sun  ^ 


then 


pert 


sun 


'pert 


sun 


sin 


RE 


«7  feet 


It  can  similarly  be  shown  that  Rpert^^„„  <  ^ 
and  lunar  gravity  can  be  ignored. 

Residual  Errors 

After  all  the  preceding  empirical  and  semi-empirical  corrections  are  made 
there  are  still  residual  errors.  These  errors  are  from  three  primary 
sources: 

a.  Inaccuracies  in  the  Approximating  Functions 

Inevitable  approximations  exist  in  any  empirical  technique. 
This  is  especially  true  for  the  functions  given  in  this 
section,  since  they  are  primarily  a  feasibility  demon¬ 
stration  and  could  certainly  be  improved. 

b.  Neglected  Effects 

All  terms  except  the  J  term  in  the  potential  expansion  are 
examples . 

c .  Neglect  of  the  Inter-Dependence  of  the  Various  Effects 

Trajectory  perturbations  due  to  oblateness  obviously  have 
an  effect  on  re-entry  perturbations  and  vice  versa. 

Likewise  trajectory  dispersions  due  to  non-nominal  missile 
parameters  affect  the  oblateness  correction  (but  not  re- 
^try  since  r,_  is  a  function  of  V„).  It  would  be  possible 
to  correct  for  these  effects  empiHcally.  although  they  are 
quite  small. 
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A  r6pr« tentative  plot  of  these  residual  errors  as  determined  from  simulation 
is  given  in  Figure  3.  It  would  of  course  be  possible  to  tabulate  these  resid- 
ual  errors  and  utiliae  then  to  obtain  extreme  accuracy  if  a  small  amount  of 
prelaunch  computation  is  acceptable. 


APPENDIX  A 


GUIDANCE  EQUATIONS 


Equations  which  can  be  used  to  generate  thrust  termination  and  steering 
commands  are  discussed  in  detail  in  Reference  1.  The  derivation  is  lengthy 
and  will  not  be  given  here.  These  equations  are  written  without  reference  to 
any  nominal  missile,  which  means  that  dispersions  due  to  missile  parameter 
variations  vanish  except  for  secondary  variations  in  the  non-Keplerian  effects. 
The  vehicle  is  steered  to  obtain  V  =  Vy  =  VyD  and  r=  at  burnout, 

where  the  desired  quantities  are  defined  from  Sections  3  and  4. 

An  illustrative  example  of  the  essential  elements  of  a  guidance  equation 
mechanization  is  given  in  the  following  figures.  The  symbol  t  represents 
the  duration  of  the  computing  cycle.  Any  value  of  t  between  0.  5  and  1.  5 
seconds  will  work  equally  well  for  practical  purposes. 
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Figar*  A-1 .  Pralaoach  Coo^vtatioM 
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Figure  A -2.  Coordinate  Transformation 
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Fig.  A-3 


Figure  A *4.  Perturbation  Corrections 
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